A general coupled nonlinear oscillator model for event-related (de) synchronization 
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Changes in the level of synchronization and desynchronization in coupled oscillator systems due to 
an external stimulus is called event related synchronization or desynchronization (ERS/ERD). Such 
changes occur in real life systems where the collective activity of the entities of a coupled system 
is affected by some external influence. In order to understand the role played by the external 
influence in the occurrence of ERD and ERS, we study a system of coupled nonlinear oscillators 
in the presence of an external stimulus signal. We find that the phenomena of ERS and ERD are 
generic and occur in all types of coupled oscillator systems. We also find that the same external 
stimulus signal can cause ERS and ERD depending upon the strength of the signal. We identify 
the stability of the ERS and ERD states and also find analytical and numerical boundaries between 
the different synchronization regimes involved in the occurrence of ERD and ERS. 

PACS numbers: 



I. INTRODUCTION 

Synchronization is an ubiquitous natural phenomenon 
that occurs widely in real systems including those in 
physics chemistry 0-0], biology 0-0] and nano- 
technology [1, This phenomenon is an active topic 
of current research and is being extensively studied [lOr - 
Il4j ]. Nevertheless, synchronization is not always a desir- 
able phenomenon. In some cases it is desirable while in 
some other times it is undesirable and hence a mecha- 
nism to desynchronize becomes necessary for normal be- 
havior. For example, synchronization is desirable in the 
cases of lasers and Josephson junction arrays [HI, Il6j . 
coupled spin torque nano-oscillators where coherent mi- 
crowave power is needed (l7l Il8j , and in the brain when 
synchronization of neuronal oscillations facilitate cogni- 
tion via temporal coding of information [l9l - [2lj | . On the 
other hand, synchronization is undesirable when pedes- 
trians walk on the Millennium Bridge [22| and when mass 
synchronization of neuronal oscillators occurs at a partic- 
ular frequency band resulting in pathologies like trauma, 
Parkinson's tremor and so on [23| . Apart from these, 
synchronization or desynchronization of neuronal oscil- 
lations in the brain are found to facilitate task selection 
and performance. When such a synchronization or desyn- 
chronization occurs in the brain or other physical systems 
mentioned above due to a particular task or an event or 
an influence of an event then it is called an event-related 
synchronization (ERS) or desynchronization (ERD). 

In this direction in order to understand the phenom- 
ena of ERD and ERS, we study the synchronization and 
desynchronization dynamics of a system of coupled non- 
linear oscillators due to the effect of the strength of the 
external stimulus which represents some external influ- 
ence. In particular, we demonstrate the occurrence of 
ERD and ERS due to the change in the strength of 
the external stimulus. We also find the occurrence of 
ERD/ERS in periodic and chaotic systems and also with 
different forms of coupling. We therefore find that the 
occurrence of ERD and ERS due to external influence is 



generic and is not limited to the type of the system or 
the coupling. 

The plan of the article is as follows. In the following 
section we present a general model of coupled nonlinear 
oscillators in the presence of an external influence to ex- 
plain ERD/ERS. We discuss the case of the limit cycle 
Stuart-Landau (SL) oscillators in Sec. IIIII as an illustra- 
tion of our model. We quantify the strength of ERD /ERS 
by defining the intensity I and identify the analytical sta- 
bility boundaries of the ERD and ERS states. Further we 
find out the stability of the ERD and ERS states and also 
identify the analytical and numerical bifurcation bound- 
aries between the different synchronization regimes. We 
also demonstrate the occurrence of ERD and ERS in the 
presence of other forms of couplings like weighted cou- 
pling and time delayed coupling. In Sec. IIVI we show 
that ERD/ERS phenomena also occur in chaotic systems 
by illustrating the case of coupled Rossler oscillators. We 
point out potential applications of the model and the re- 
sults in Sec. [V] Finally we present our conclusions in 
Sec. ED 



II. THE GENERAL MODEL FOR ERD/ERS 

Let us consider a system of coup led limit cycle non- 
linear oscillators 0, [LL, [H subject to external 
stimulus. The external stimulus represents a dynamic 
signal from outside the system (brain or or other physi- 
cal systems such as spin torque nano oscillators) or from a 
distant region inside the system and typically represents 
the influence of the event (and not necessarily the ac- 
tual event itself). With this in consideration, the model 
equations shall be cast into the following form, 

A N d 
X, = F(X j ,e j ) + -J2(*k-X j ) + BY, (•=-) 



k=l 



G(Y,e e ), 



(1) 
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where j = 1, 2, . . . , N. Here F(Xj, ej) represents the non- 
linear limit cycle behavior of the jth uncoupled oscillator 
and 6j is the corresponding system parameter. A is the 
coupling strength between the oscillators in the system 
and B represents the coupling strength between the os- 
cillators in the system and the external stimulus. Our 
findings on the occurrence of ERD in this model has been 
briefly reported in |28[ . 

From a detailed analysis of generic models of the form 
(TTJ) , we have found that the external stimulus B typically 
affects synchronization in the system of coupled oscilla- 
tors. If the system of coupled oscillators is completely 
synchronized, that is, X 3 = X, j = 1, 2, . . . , N (in this 
state all the oscillators behave as one) due to the strength 
of the coupling A (when B = 0), then a sufficient strength 
of B causes desynchronization (ERD, since external stim- 
ulus causes desynchronization) in the system. Typically 
ERD occurs due to a small group of oscillators that sepa- 
rate themselves from the larger synchronized group. The 
separated small group either remains desynchronized or 
synchronizes itself to a frequency different from that of 
the larger synchronized group. Both of these cases will 
be demonstrated in the following Sections. 

The above behavior is generic and can also oc- 
cur in chaotic systems leading to chaotic synchroniza- 
tion/desynchronization. Further, this behavior occurs ir- 
respective of the form of coupling (weighted, nonlocal, 
etc.) and also in delay coupled systems. In all the cases, 
when the system is synchronized due to the strength of 
the coupling and when the strength of the external stim- 
ulus is sufficient, ERD/ERS occurs. In the following we 
present our theoretical findings and discuss the qualita- 
tive connection between our findings and experimental 
observations. 



III. ERD/ERS IN A SYSTEM OF 
STUART-LANDAU OSCILLATORS 

A well-known and well-studied limit cycle nonlinear os- 
cillator is the Stuart-Landau (SL) oscillator [4]. A system 
of coupled SL oscillators is described by the following set 
of coupled complex first order nonlinear ordinary differ- 
ential equations (ODEs), 



A 



N 



(a + icjj - (I + ic)\Zj\ 2 )zj + — ^2(z k - Zj) + B 



k=l 



z e = (a e + iuj e - (I + ic e )\z e \ 2 )z e 



(2) 



Here c is the nonisochronicity parameter, a is the Hopf 
bifurcation parameter and T^j = (z^ — Zj) is the coupling 



function. Zi 



+ iyj is the complex amplitude of the 



jth oscillator with natural frequency uij . z e is another SL 
oscillator with natural frequency u e and nonisochronicity 
parameter c e . Without coupling (A, B — 0), the dynam- 
ical equation of an individual Stuart-Landau oscillator 
becomes 





MB 



r = (a — r 2 )r, 



w — cr 



(3) 



FIG. 1: Occurrence of ERD and ERS in a system of 1000 
coupled Stuart-Landau oscillators for A = 1.1, c = 1.5, c e = 
2.0, a = a e = 1.0, cj e = 1.5,7 = 0.005 and u) — 0.5 in equation 
([2]). Left column shows the time evolution of the real part 
of the state vectors, Xi, and the right column depicts the 
corresponding phase portraits on the (xi,yi) plane. Here (top) 
B = 0.01, (middle) B = 0.08 and (bottom) B = 0.35. Note 
%i — %i ~\~ yt- 



where z = re 19 . The fixed points of the radial equation 
are ro = and y/a. From a linear stability analysis the 
eigenvalue is found as A = (a— 3rg). For a < the system 
has a stable fixed point z — which is a linearly stable 
state corresponding to amplitude death. For a > 0, the 
fixed point z — becomes linearly unstable and a limit 
cycle oscillation is established with radius r — yfa and 
phase 9 — (u>—ac)t+6o, where #o is a constant. For more 
detailed studies about synchronization in populations of 
SL oscillators one may refer to [HI, HH . 

The occurrence of ERD and ERS in a system of 1000 
coupled Stuart-Landau oscillators is numerically demon- 
strated in Figure [T] for a Lorentzian distribution of nat- 
ural frequencies g(u>) = ^[y 2 + (uj — cj) 2 ] -1 , where 7 is 
the half width at half maximum and u> is the central fre- 
quency. For the given set of system parameters (A = 1.1, 
c = 1.5, c e = 2, u e = 1.5, 7 = 0.005, ui = 0.5), ERD oc- 
curs when moving from the top to the mid panel and 
ERS occurs when moving from the mid to the bottom 
panel, upon increasing B. The time evolution of the real 
part of the state vectors, and the corresponding phase 
portraits are plotted in the left and the right columns, 
respectively. When B = 0.01 all the oscillators are syn- 
chronized in phase (while there is slight desynchroniza- 
tion in the amplitude) in the top panel. On increasing 
B from 0.01 to 0.08, in the mid-panel, ERD occurs when 
some oscillators separate themselves from the synchro- 
nized group and are synchronized to a different frequency. 
Even though the separated group of oscillators are syn- 
chronized among themselves, we refer to this state as 
being desynchronized because from the point of view of 
the entire system, this is desynchronization. Further due 
to this separation, the intensity or the strength of syn- 
chronization (discussed in detail in the following subsec- 
tion) of the original synchronized group (shown in the top 
panel) is reduced. Now, in this state, when we increase 
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FIG. 2: The case of quasi-periodic synchronization in the 
ERD state in a system of 1000 coupled Stuart-Landau oscil- 
lators. Here ui e = 1.0; all the other parameter values are the 
same as in Figure [1] 



FIG. 3: Change in the intensity I for the SL system (J2j f° r 
varying stimulus strength B for different values of the cou- 
pling strength A. All other parameters are the same as in 
Figure [2] In the direction of increasing B, downward arrow 
(that is, a decrease in 7) denotes ERD and the upward arrow 
(that is, an increase in /) denotes the occurrence of ERS. 



B from 0.08 to 0.35 ERS occurs (with reference to the 
mid-panel) when the separated group of oscillators again 
synchronize with the other oscillators. The separated 
group of oscillators need not always be synchronized but 
can be either desynchronized or quasi-periodically syn- 
chronized. The case of quasi-periodic synchronization is 
shown in Figure [2j Here ERD occurs in the mid-panel for 
the same values of (system and control) parameters as in 
Figure [TJ except that now u e = 1.0, where the separated 
group is quasi-periodically synchronized. 



A. Measure of ERD/ERS 

The strength or intensity of synchronization can be 
measured by the number of oscillators that are oscillat- 
ing in synchrony. The more the number of oscillators 
that are oscillating in synchrony the more will be the 
strength of synchronization. Hence an increase in the 
strength or intensity of synchronization denotes the oc- 
currence of ERS and a decrease in the intensity denotes 
the occurrence of ERD. With this reasoning we quantify 
the intensity of synchronization using the phases of the 
oscillators by denning the quantity [J HH 



I=< 



i 

>= T 



1 N 



3 = 1 



dt. 



(4) 



where 9j — taxT l {yj/xj) is the phase of the jth oscil- 
lator. Here the bar represents average over all oscilla- 
tors in the population and the angle brackets represent 
time average. In the state of complete desynchronization 
1 = and for complete synchronization 1=1 (while we 
neglect very little desynchronization in the amplitude). 
For partial synchronization / takes a value between 
and 1; the more oscillators are oscillating in synchrony 
the higher will be the value of /. We use phase desynchro- 
nization/synchronization to characterize the occurrence 
of ERD/ERS because, upon increasing B while the os- 



cillators are in complete synchronization, very little am- 
plitude desynchronization occurs first and then for an 
increase in B desynchronization occurs in the phase and 
the desynchronization in the amplitude is also increased 
(ERD) (see for example Figures Q] and [2]) . Further in- 
crease in B brings back synchronization in the phase with 
slight desynchronization in the amplitude (ERS). There- 
fore monitoring the intensity of phase synchronization 
will facilitate monitoring the occurrence of ERD/ERS. 

Figure [3J shows that there exists a critical strength of 
the external stimulus for ERD/ERS to occur. Upon vary- 
ing the stimulus strength B, ERD occurs when there is 
a decrease in the intensity at a sufficient B (as denoted 
by a downward arrow in Figure [3]) and when there is an 
increase in the intensity for sufficient B, ERS occurs (as 
denoted by an upward arrow in Figure [3j. It may also be 
noted that the ERD/ERS occurs for a finite window of 
the stimulus strength which depends upon the coupling 
strength A. When we start with a state where the sys- 
tem is completely phase synchronized due to the coupling 
strength A, the ERD/ERS occurrence window decreases 
with stronger A. This means that for very strong cou- 
pling A there is a very less chance that ERD/ERS occurs. 



B. Stability of the ERD/ERS state 

Our numerical simulations show that when ERD oc- 
curs for a sufficient strength of the stimulus, the syn- 
chronized group splits into one major synchronized group 
and a small separated group (synchronized or desyn- 
chronized) comprising of few oscillators (say N^) com- 
pared to the size of the major synchronized group (say 
Ni ) . Therefore the size ratio of this two cluster state is 
r : 1 — r where r = Ni/N and N = Ni + N%, and in 
the desynchronized state r >> (~ 1). For the case of 



SL oscillators © with an external force 



(so 



that c e = and a e = 1), the dynamics of the large 



4 



group Zj 



u, J 



1,2,...N\ and the small group 



Zj = v, j = Ni + l,Nx + 2, . . . N can be written as 

u fa ((a + iu) - (1 + ic)\u\ 2 )u + Be l ^\ 

v fa ((a- A + iuj) - (1 + ic)\v\ 2 )v + Au + Be iUet .(5) 

If we assume that the major synchronized group is com- 
pletely synchronized with the external force z e = e lWet , 
while the small group is not, we get u> = uj e + c and 
a = 1 — B. From this reasoning, the dynamics of the 
small group can be written in the form 

w = ((a - A + ic) - (1 + ic)\w\ 2 )w + A + B, (6) 

where w = ve~' lu " t . The occurrence of ERD/ERS can be 
understood from the fixed points of this equation [28j . 
Due to the cubic nature of equation ([6]) it has three fixed 
points which determine the stability of the ERD/ERS 
state. One of the fixed points, namely w\ = 1, is stable. 
The eigenvalues of the linearized version of (jSJ) are \\ t 2 = 
— (A + B),—(2 + A + B). This corresponds to a complete 
synchronization (synchronization of all the oscillators) of 
the population, that is u = v [HJ . The other two fixed 
points (say W2 and Ws) determine the stability of the 
desynchronized state. The fixed points W2 and W3 exist 
for 



B < B T = (1 + c 2 )(y/l + c 2 - l)/(2c 2 ) - A, (7) 

where Bj is a saddle-node bifurcation point and the fixed 
points are given by 



1 



2(1 + c 2 ) 



1 + c 2 - 2(A + B) 



±y/(l + c 2 - 2{A + B)) 2 - 4(1 + c 2 ){A + B)m 



W2,3 = 



Performing a linear stability analysis on equation ([6]) we 
find that the determinant of the Jacobian matrix for W2 
and W3 is 



det(J W23 ) 



T2{A + B)(A + B + 2)A 
{1 + c 2 + 2{A + B))±A ' 



(9) 



where A = ((1 + c 2 + 2(A + B)) 2 -4(1 + c 2 ){A + B){A + 
B + 2)) 2. It is easily checked that the determinant (O 
is always negative (positive) for u> 2 (w^). Therefore the 
fixed point W2 turns out to be a saddle. The other fixed 
point W3 is either an unstable node or focus for |c| < 1. 
When \c\ > 1 and B > Bjj, the fixed point W3 is either a 
stable node or focus; here Bu is a Hopf bifurcation point 
given by 



B u = (1 + c 2 )/(v/4 + (1 + c 2 ) 2 + 2) -A, (10) 
which is determined from the condition tr( J W3 ) = 0. Here 
(l-c 2 )(A + B) + A 



tr(J W3 ) = 2- 



(1 + c 2 ) 



(11) 
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FIG. 4: B— |c| phase diagram. Bi (solid line) and Bu (dashed 
line) are the analytically obtained saddle-node and Hopf bifur- 
cation boundaries respectively. The boundary Bin (dotted 
line) is obtained numerically by solving the evolution equation 
for w given in the text. 



Further for B < Bm only synchronized solutions are 
stable and desynchronized solutions do not exist. For 
given values of parameters the bifurcation boundaries 
are plotted in Figure |4] In the regions I and IV all 
the oscillators are synchronized, while in the region II 
there is frequency desynchronization (mid-panel of Fig- 
ure [2] is an example since the small group oscillates quasi- 
periodically while the major group oscillates periodically) 
between the two clusters and in the region III there is am- 
plitude desynchronization between the two clusters. At 
the point £?///, both the saddle and the Hopf bifurca- 
tion points merge and disappear and this point is called 
a saddle-connection point. It is obtained numerically by 
solving equation ([6]) and is shown as a dotted line in Fig- 
ure |4] Thus desynchronized solutions exist in the region 
between Bi and Bm. This is in agreement with our 
numerical observations. A movement in the parameter 
space in the direction / — > II or IV — > III represents 
ERD and a movement in the opposite direction repre- 
sents ERS. 



C. ERD/ERS in systems with other forms of 
coupling 

In order to demonstrate further the validity of the 
model, we first consider a system of delay coupled SL 
oscillators described by 



Zj = (a + iujj - (1 + ic)\z 3 \ 2 ) Z j 



A 



N 



-J2(z k (t-T)- Zj )+Bz e , (12) 



fc=i 



where r is the time delay parameter. We simulate (|12[) 
with t — OA and find the occurrence of ERD/ERS as de- 
picted in Figure|5](a)-(c) for A = 0.035, c = 2.5, c e = 3.0, 
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FIG. 5: The time evolution of the state vectors Xi of a system 
of 100 delay coupled SL oscillators showing the occurrence 
of ERD/ERS in panels (a) - (c). Here r = 0.1, A = 0.035, 
c = 2.5, c e = 3.0, a = a e = 1.0, u; e = 1.5,7 = 0.0005 and 
lj = 1.5 and B = 0.0 for (a), B = 0.05 for (b) and B = 0.4 for 
(c). Occurrence of ERD/ERS in a system of 100 coupled SL 
oscillators with random weighted coupling with weight factor 
3.5 is shown in panels (d)-(f). Here c = 2.5, c e = 1.0, oj e = 
0.5,7 = 0.005 and lj = 1.5 and B = 0.0 for (d), B = 0.1 for 
(e) and B = 1.1 for (f). 



cj e = 1.5,7 = 0.0005 and Q = 1.5. In panel (a) the oscil- 
lators in the system are synchronized due to the strength 
of coupling A when the external stimulus is absent. In 
panel (b) when the external stimulus is switched on with 
a strength B = 0.05 ERD occurs in the system where a 
few of the oscillators separate themselves from the major 
synchronized group. When the strength of the external 
stimulus is further increased to B = 0.4 in panel (c) ERS 
occurs in the system and all the separated oscillators get 
back to synchronization with the major group. Secondly, 
ERD /ERS also occurs in a system with random weighted 
coupling described by 



zj = (a + iuij - (1 + ic)\zj\ 2 )zj 

N 

+ ^A jk {z k - Zj )+Bz e , (13) 
k—l 



where A\j is the random coupling matrix. Figure [5] 
(d)-(f) demonstrates the occurrence of ERD (panel (e)) 
and ERS (panel (f)) in this system due to the stimulus 
strength. Thus we find that the form of the coupling is 
immaterial in causing the occurrence of ERD/ERS. For 
the non-delayed case, if we reduce our model to a phase 
model the bifurcations match with those presented in 
[35H37j where the authors have studied a system of cou- 
pled phase oscillators in the presence of external forcing. 



IV. ERD/ERS IN A SYSTEM OF COUPLED 
CHAOTIC OSCILLATORS 

The occurrence of ERD /ERS is generic and not limited 
to limit cycle oscillator systems alone. For example let 
us consider a system of coupled Rossler oscillators in the 
presence of an external field (B ^ 0) 

Xj = —Ldjl/j — Zj + Bx e , 

Vi = UjXj + ay } + — ^{y k - y 3 ), 
fc=i 

ij = b + Zj(xj — c), (14) 
where the external stimulus is described by 

x e = -uj e y e - z e . y e = u)jX e + ay e , z e — b e + z e (x e - c e ). 

Also this system is found to exhibit ERD/ERS. Here a, 
b, and c are the system parameters that determine the 
periodicity/chaoticity of the system. A Lorentzian dis- 
tribution of the natural frequencies uij are chosen for nu- 
merical simulation. The external stimulus is represented 
by x e and in this case we have introduced coupling with 
the x component of the external stimulus. However, the 
phenomenon demonstrated here can also be observed ir- 
respective of the component in which the external stimu- 
lus is coupled to. The system parameters of the coupled 
Rossler oscillators are chosen so that they are in the pe- 
riodic regimes while those of the external oscillator are 
chosen so that it is in the chaotic regime. 

An illustration of ERD and ERS in this system of 1000 
coupled oscillators is depicted in Figure|B]for a Lorentzian 
distribution of natural frequencies with ui = 1, a = 0.1, 
b = 0.5, c = 4.0 and a e = 0.1, b e = 0.1, c e = 9.0 and 
7 = 0.005. The top panel is the state of complete phase 
synchronization due to the coupling strength A = 0.7 
when B = 0.1. ERD occurs in the mid-panel for B = 1.8. 
In this state some of the oscillators separate from the syn- 
chronized group and are desynchronized. However unlike 
the case of ERD in SL oscillators (depicted in Figure [T]) 
this ERD state is not stable but breathes in time. That 
is the separated group is desynchronized at certain times 
while it is synchronized with the major group at other 
times. This is obvious from Figure [7] which is a zoom 
of the rectangular regions marked in Figure |51 Figure [7] 

(a) corresponds to the dotted rectangular region, which is 
the state where some oscillators are separated and desyn- 
chronized from the major synchronized group. Figure [7J 

(b) corresponds to the solid rectangular region, which is 
the state where all the oscillators are synchronized. We 
have confirmed that this is not a transient since we have 
plotted the data after allowing a sufficiently large number 
of iterations and the breathing phenomenon is found to 
occur periodically. Thus ERD in this state is not stable 
but time varying. On further increase in the strength of 
the stimulus to B = 2.3, ERS occurs in the system when 
all the oscillators are synchronized. 
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FIG. 6: Occurrence of ERD and ERS in a system of 1000 
coupled Rossler oscillators. Here A = 0.7, a = 0.1, b — 0.5, 
c = 4.0 and a e = 0.1, b e = 0.1, c e = 9.0, 7 = 0.005, Q = 1.0 
and (top) B = 0.1, (mid) B = 1.8 and (bottom) B = 2.3. Left 
column shows the time evolution of the state vectors Xi and 
the right column depicts the corresponding phase portraits 
on the (xi,yi) plane. A blow up of the dotted and solid rect- 
angular regions marked in the mid-panel is shown in Figure 
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FIG. 8: Occurrence of ERD in a system of 1000 coupled 
Rossler oscillators that are in the chaotic regime. In panels (a) 
and (b) the external oscillator is also chaotic with a e = 0.2, 
b e — 0.2, c e = 5.7. In panels (c) and (d) the external oscilla- 
tor is periodic with a e = 0.1, b e = 0.5, c e = 4.0. ERD occurs 
in panels (b) and (d). Here A = 0.7, a = 0.2, b = 0.2, c = 5.7, 
7 = 0.005, Q = 1.0 and B = 0.01 for (a), B = 0.15 for (b), 
B = 0.03 for (c) and B = 0.21 for (d). 
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FIG. 7: A blow up of the dotted and the solid rectangular 
regions marked in Figure [6] is shown in panels (a) and (b), re- 
spectively. In panel (a) the separated group is desynchronizcd 
while in (b) it is synchronized with the major group. 



The same stimulus can result in ERS and ERD de- 
pending upon the strength of the stimulus. Fig. [H] shows 
this phenomenon where we choose the initial state of the 
system to be in a partially synchronized state (since we 
want to demonstrate both ERS and ERD). Panel (a) rep- 
resents the initial state before the application of the stim- 
ulus. When the stimulus is applied, depending upon the 
strength of the stimulus B ERD occurs in panel (b) while 
ERS occurs in panel (c). The change in intensity / for 
different values of stimulus strength B is shown in panel 
(d). The reference line corresponds to the initial partial 
synchronized state and an intensity / above and below 
the line denotes the occurrence of ERS and ERD, respec- 
tively. Thus we find that for a given external stimulus 
signal, both ERD and ERS can occur depending upon 
the strength of the stimulus. 



V. APPLICATIONS 



It may be noted that for illustrative purpose we have 
shown here the occurrence of ERD/ERS in a system of 
coupled Rossler oscillators that are in periodic regime 
stimulated by a chaotic Rossler oscillator. However, the 
same can be observed even if the external oscillator is 
in the periodic regime and coupled oscillators are in the 
chaotic regime and also for cases where the coupled os- 
cillators and the external oscillator are both in chaotic 
or periodic regimes. For illustration we have shown the 
occurrence of ERD in a system of coupled chaotic Rossler 
oscillators with a chaotic stimulus in panels (a) and (b) 
and with a periodic stimulus in panels (c) and (d) of 
Figure |U 



In this paper we have discussed the occurrence of 
event-related synchronization and desynchronization in 
a system of coupled nonlinear oscillators in the presence 
of an external stimulus. This type of model and the re- 
sults discussed here can be applied to various real world 
systems. We point out a few applications below. 

Synchronization of neurons at a particular frequency 
band may facilitate the performance of a particular task 
while desynchronization at the same frequency band may 
leave it unattended. The vice verse can also happen. 
Brain oscillations at different frequency bands are found 
to be one of the most crucial mechanisms that control 
higher level information processing, motor functioning 
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FIG. 9: Numerical illustration of the occurrence of ERD and 
ERS in a system of 1000 SL oscillators for different stimulus 
strengths. Here A = 0.5, c = 1.5, c e =2.0, uj e = 1.0, 7 = 0.05, 
lj = 2.5 and B = 0.0 for (a), B = 1.5 for (b) and B = 3.0 for 
(c). (a), (b) and (c) depict the phase portraits in the (xi,yi) 
plane, (d) shows the change in the intensity / for different 
values of B, representing ERD and ERS. 



ERD/ERS phenomena in the brain, can also explain syn- 
chronization and/or desynchronization due to external 
stimulus in other systems; for instance, polariton conden- 
sates in semiconductor micro-cavities that interact both 
among themselves and with the reservoir and quantum 
coherence of condensates [29|, [3(| can also be represented 
by this model. Studies on these systems will be discussed 
separately. 

Even more generally, the model can account for syn- 
chronization/desynchronization due to external stimulus, 
while the system is already synchronized due to coupling. 
Possible applications in this direction include the estab- 
lishment of synchronization or desynchronization due to 
external stimulus in neuronal networks, Bose-Einstein 
condensates, Josephson junction arrays and lasers and 
so on. 

Therefore it will also be of interest to investigate the 
occurrence of ERD/ERS in systems of diffusively coupled 
sub-populations. Such a model will facilitate the under- 
standing of the simultaneous occurrence of ERD/ERS at 
multiple frequency bands in the brain and in systems of 
coupled STNOs in multiple columns. 



and even large scale integration of information across var- 
ious regions of the brain. 

It is therefore highly important to understand the dy- 
namics of event related neuronal oscillations in relation to 
specific tasks and/or pathologies. In particular, in patho- 
logical conditions like tremors, it is of great importance to 
understand the occurrence of synchronization or desyn- 
chronization of neuronal oscillations at specific frequency 
bands in order to be able to control the synchronization 
of neuronal oscillations and hence the pathological states. 

The occurrence of synchronization or desynchroniza- 
tion due to an event is not limited to brain and there are 
many other systems that exhibit this phenomenon. For 
instance, an open problem in the field of nano-technology 
is that the microwave power emitted by a single spin 
torque nano-oscillator (STNO) is very small, of the or- 
der of nano-watts. In order to increase the power one 
needs to find ways to synchronize a group of coupled 
nano-oscillators to obtain coherent power. The model 
discussed in this paper can be used to find ways to avoid 
desynchronization in such systems and can also explain 
the occurrence of synchronization and/or desynchroniza- 
tion caused by a stimulated microwave current (l7| . 

We also note here that this model, apart from serv- 
ing as a qualitative mathematical representative for 



VI. CONCLUSION 

In this article we have proposed a general model of 
coupled nonlinear oscillators which can represent inter- 
esting physical or biological systems such as a system of 
coupled spin torque nano-oscillators or neurons, by gen- 
eralizing and building upon our previous Letter [28| . We 
subject the system to an external field, which represents 
an event or an external stimulus. The model results pro- 
vide an understanding of the occurrence of ERD/ERS 
due to different types of stimuli. We have found that the 
occurrence of ERD/ERS is generic to all types of cou- 
pled oscillator systems including chaotic systems. Hence 
we hope that the model results can be applied to quan- 
titative description and prediction of the occurrence of 
synchronization and desynchronization due to external 
stimulus in real world systems including Bose-Einstein 
condensates, Josephson junction arrays and lasers and 
so on. 
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